Structure of the breakpoint region in CVC of the intrinsic Josephson junctions 
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A fine structure of the breakpoint region in the current-voltage characteristics of the coupled 
intrinsic Josephson junctions in the layered superconductors is found. We establish a correspondence 
between the features in the current-voltage characteristics and the character of the charge oscillations 
in superconducting layers in the stack and explain the origin of the breakpoint region structure. 



A system of Josephson junctions attracts a great inter- 
est from both scientific and practical sides. To describe 
it, a system of coupled nonlinear equations is used, which 
might be solved numerically only. Recently observed^ 
powerful enough coherent radiation from the stack of the 
intrinsic Josephson junctions (UJ) in layered supercon- 
ductor BSCCO opens new perspectives for different ap- 
plications. 

In Refsi^i 3 -'^ we studied the multiple branch structure 
of the current-voltage characteristics (CVC) of IJJ and 
showed that the branches have the breakpoint (BP) and 
a breakpoint region (BPR) before transition to the an- 
other branch. The BP current is determined by the cre- 
ation of the longitudinal plasma waves (LPW) with a 
definite wave number k, which depends on the coupling 
parameter a, the dissipation parameter /3, the number 
of junctions in the stack, and the boundary conditions. 
We generalized the McCumber-Stewards dependence of 
the return current for the case of IJJ in the HTSC and 
the oscillation of the breakpoint current as a function 
of parameters a and (3 was found. Based on the idea 
of the parametric resonance in the stack of IJJ, we ex- 
plained this oscillation as a result of the creation of lon- 
gitudinal plasma waves with different wave numbers. A 
good qualitative agreement of the results of modeling of 
the a/3-dependence of the BP current with the results 
of simulation was obtained. We demonstrated that the 
a/3-dependence of the BP current is an instrument to de- 
termine the mode of LPW created at the breakpoint in 
the stacks with a different number of junctions. So, the 
breakpoint in the CVC was explained, but the nature 
and features of the breakpoint region were not clear at 
that time. 

In this paper we study theoretically the phase dynam- 
ics of a stack of IJJ in the high-T c superconductors. The 
CVC of IJJ are numerically calculated in the framework 
of capacitively coupled Josephson junctions model with 
diffusion current From the results of the numerical 
simulations we predict a fine structure of the breakpoint 
region in the CVC. A correspondence between the fea- 
tures in the CVC and the character of the charge oscil- 
lations on different superconducting layers is established 
and the origin of the breakpoint region structure is ex- 
plained. In the BPR the plasma mode is a stable solution 
of the system and this fact might be used in some applica- 
tions, particularly, in high frequency devices such as THz 



oscillators and mixers. We consider that the breakpoint 
phenomenon might be used to understand the mechanism 
of the powerful coherent radiation from the stack of the 
IJJ in BSCCO. 

To find the CVC of the stack of IJJ, we solve a sys- 
tem of dynamical equations for the gauge- invariant phase 
differences tpi(r) = 0j+i(t) - 0j(r) - ^ dzA z (z,r) 
between superconducting layers (S-layers) for the stacks 
with a different number of IJJ in the framework of the 
capacitively coupled Josephson junctions model with dif- 
fusion current (CCJJ+DC model)^^, where Oi is the 
phase of the order parameter in the S-layer I, A z is the 
vector potential in the barrier. We use a dimension- 
less time r = tw p , where oj p is the plasma frequency 
Lu p = yj2eI c /hC , I c is the critical current and C is 
the capacitance. The system of equations has a form 
d 2 ipi/dr 2 = J2i> Aw[I — sin</j;/ — fid<pv /dr] with matrix 
A given in Ref.— for periodic and nonperiodic boundary 
conditions (BC). In our simulations we measure the volt- 
age in units of Vq — hw p /{2e) and the current in units 
of the I c - The CVC and time dependence of the charge 
oscillation in the S-layers are simulated at a = 1 and 
j3 = 0.2 and periodic BC. The details concerning the 
model and numerical procedure are presented in Refsi^£ 
In the inset to Fig. we show the result of simulation 
of the CVC for a stack with 9 IJJ. It demonstrates an 
equidistant multiple branch structure in the hysteresis 
region and the breakpoint location. The transition from 
the state where the outermost branch with all junctions 
is in a rotating state (R-state) 7 to a state with some 
junctions in an oscillating state (O-state ) happens in a 
resonance region, where the LPW with a definite wave 
number is created. In the stack with 11 IJJ with men- 
tioned above values of a = 1 and (3 — 0.2, the LPW with 
k = 107r/lld is created, where d is the lattice period 
along the c-axisi^ In Fig. QJ, we can see that the CVC 
exhibits a fine structure in the BPR. The BPR has dif- 
ferent parts with a different dependence of the voltage on 
the bias current. For a given stack we can see clear five 
parts in the CVC separated by some points T, Y, C\,Ci, 
where the character of the CVC is changed. The letters B 
and J show the breakpoint and jumping point to the an- 
other branch. The points S on the CVC will be discussed 
later. 

To find the origin of this structure and explain its 
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FIG. 1: (Color online) (a) The fine structure of the BPR in 
the outermost branch in the CVC of the stack with 11 IJJ. 
The insert shows a total branch structure in the CVC of the 
stack with 9 IJJ and the BP location, (b) Charge oscillation 
in the first layer in the beginning of B-S part of the BPR. The 
inset shows the oscillations in the layers 1 and 2. 

features, we study time dependence of the charge in 
the superconducting layers. Using Maxwell equation 
div(eeoE) = p, where e and Eq are relative dielec- 
tric and electric constants, we express the charge den- 
sity Qi (we call it just charge) in the S-layer I by the 
voltages Vi and Vi+i in the neighbor insulating layers 
Ql = Qoa(Vi+i - Vi), where Q = eeoV /r 2 D , and r D 
is Debay screening length. Solution of the system of 
dynamical equations for the gauge-invariant phase dif- 
ferences between S-layers gives us the voltages Vi in all 
junctions in the stack, and it allows us to investigate the 
time dependence of the charge on each S-layer. 

The "time dependence" actually consists of time and 
bias current variation. We solve the system of dynami- 
cal equations for phase differences at fixed value of bias 
current / in some time interval (0,T m ) of dimensionless 
time r = tuj p with the time step St, where t is a real 
time. This interval is used for time averaging procedure. 
Then we change the bias current by SI, and repeat the 
same procedure for the current I + SI in new time inter- 
val (T m , 2T m ). The values of the phase and derivative of 
phase in the end of the first time interval are used as ini- 
tial conditions for second time interval and so on. In our 
simulations we put T m = 250, St = 0.05, SI = 0.0001 and 



total recorded time was calculated as r + T m (I — I) /SI, 
where Iq is an initial value of the bias current for time 
dependence recording. 
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FIG. 2: (Color online) (a) Profile of time dependence of the 
charge in the S-layers in the stack with 9 IJJ: (a) In the first 
layer, (b) In the third layer, (c) In the fourth layer. Red line 
shows the CVC of this stack in the BPR. 

In the Fig. [2] the profile of time dependence of the 
charge in the S-layers in the stack with 9 IJJ is combined 
with the CVC of the outermost branch at a = 1, (3 = 
0.2 and periodic BC. We show oscillations in the first 
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(a), third (b) and fourth (c) layers only, because in the 
other layers its character close to the behavior in these 
layers. To avoid large data files in the procedure of time 
dependence recording and to present the CVC and time 
dependence with the same precision, we have limited the 
averaging time interval by T m = 250. It leads to more 
smeared CVC in this figure, than the one shown in Fig.[T]a, 
(T m = 25000), but the main features of the CVC are 
preserved. Moreover, such an averaging also allows us to 
find one more characteristic point S within the BPR. 
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FIG. 3: (Color online) Superposition of the three CVC with 
the amplitudes of noise in the current 10" 6 , 1CT 8 and 10~ 13 . 



We have investigated the influence of the noise in bias 
current on the CVC, which is generated arbitrary in each 
IJJs . The noise is produced by random number gener- 
ator and it's amplitude is normalized to the critical cur- 
rent value I c . Fig. [3] presents the CVC of the outermost 
branch in the vicinity of the breakpoint for three val- 
ues of noise amplitude: 10 -6 , 10~ 8 and 10~ 13 of a stack 
with 11 IJJ. We see that first four parts of the BPR from 
the point B till the point C\ are not affected by noise 
in current practically. But its strong effect has a place 
in the C\ — J region. The noise in bias current smears 
the CVC essentially here because four first parts of the 
BPR demonstrate a behavior, close to the regular one, 
but the part C\ — J shows the chaotic behavior. The 
voltage value is very sensitive to the small changes in the 
value of current here. 

As we can see in the Fig. [2 the features of the CVC 
are in correlation with the features of time dependence 
of the charge on the S-layer. This fact can be explained 
by the idea of the parametric resonance at point B. As 
was shown in Ref£, the system of equations for CCJJ 
has a solution corresponding to the LPW propagating 
along the c axis. A frequency of the LPW at / = and 
f3 = is u>Lpw(k) = ujpy/l + 2a(l — cos fed), where k 
is wave vector of the LPW. At point B the Josephson 
oscillations excite the LPW by their periodical actions. 
The frequency of Josephson oscillations is determined by 
the voltage value in the junction, so at luj = 2u>lpw the 



parametric resonance is realized and the LPW is created. 

The sharp increase of the oscillation amplitude in the 
first layer of the stack with 9 IJJ in the beginning of the 
B-S range is demonstrated in Fig. [TJa, where the LPW 
with wave vector k = 8ir/9d is created. This mode is close 
to the 7r-mode (k — n/d). In the 7r-mode charges on the 
nearest lattice points oscillate in the opposite directions 
and the inset shows by numbers 1 and 2 that the charges 
oscillate in the layers 1 and 2 in opposite directions too. 
But the amplitude of oscillations in the layer 2 is smaller 
than in the layer 1, which demonstrate the difference of 
this mode from the tt- mode. Fast Fourier transfer- 
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FIG. 4: (Color online) Results of FFT analysis of time depen- 
dence of the charge in the third S-layer in the different parts 
of BPR 




(a) S-T, I / l c =0.569 



r (b) T -Y, I / 1 =0.565 

0.5 fji |||| III ||| ||| ||| 

-, illy mini mini hi 




2538( 

e, tCOp 



Time, tujp 
r(c) Y-C„ I/I =0.5625 




37600 37700 

Time, tCOp 

(d) C, - J, l/l =0.556 




601 00 60200 
Time, tCOp 



FIG. 5: (Color online) Charge oscillations in the third S-layer 
in the stack with 9 IJJ in different parts of the BPR, specified 
by the value of I/I c 

mation (FFT) analysis for the time dependence of the 
charge in the third S-layer presented in Fig. |4j The time 
domain samples for FFT are collected across the interval 
(0, T m ) at each step of the bias current. Fig. Hh. shows 
that before the BP at I/I c = 0.573 the Josephson fre- 



4 



quency toj = 0.4542 * 27rw p = 2.8538w p is observed only. 
In accordance with the parametric resonance, we can see 
the peak at ui = 0.2246 * 27rw p = 1A112u> p corresponding 
to the LPW frequency lolpw in the B-S region (Fig. 0b). 
We observe also the peak at u> — 0.6738 * 2ttlu p = 4.2336 
corresponding to sum of the Josephson and LPW fre- 
quencies luj + u) lpw- A sharp increase in the oscillation 
amplitude at the resonance point is observed on other lay- 
ers in the B-S region also. The oscillation in the S- T part 
is similar to those for all layers in the stack excluding the 
layer 4. Because we consider the periodic BC along c-axis 
here, the layer number is used just for distinctness. As we 
can see in Fig. [2(c), the amplitude of the charge oscilla- 
tion in the layer 4 in the B- Y region is essentially smaller 
than that in the other layers. The FFT analysis of the 
S- T part shows the same peaks as for iJ-S'part, but here 
we can see an additional peak 0.4395 * 27rw p = 2.7615u p , 
which value approximately equal to 2lolpw- As it is 
demonstrated in Fig. [5(b), the character of the charge 
oscillations in T- Y part of the BPR is different from that 
in the S- T part of the BPR. The charge oscillations are 
characterized here by the beating process. A specific fea- 
ture of these part is the absence of the charge oscillations 
in the 4-th layer. At point Y the charge oscillations in 
the layer 4 grow rapidly and are getting the same order 
as the value in the other layers. Due to the coupling 
between junctions the oscillations in the other layers are 
increasing too, for example, like in the third layer shown 
in Fig. |2(b). The Y — C\ part of the BPR (Fig. [5(c)) 
is characterized by the beating process as well, but with 
an additional frequency in compare with the T- Y region. 
The Ci — J part of the BPR presented in Fig. [5]i dis- 
play a behavior close to the chaotic one. So, the different 
character of the charge oscillation in the S-layers is in 
correlation with the features of the CVC and reflects its 
nonlinear behavior in the BPR. 

The question concerning the origin of the BPR and 
the fine structure in the CVC might be simplified if we 
take into account the results obtained for stacks with 
even number of IJJ in Ref. 3 . In that case the LPW with 
wavelength A = 2d (k = ir) is created. This wavelength of 
the LPW is commensurable with the length of the stack 
in the c-direction and the BPR in the CVC does not 
appear. For such stacks with even number of junctions 
the resonance region is observed only, which corresponds 



to the continuous increase of the amplitude of the charge 
oscillations and followed by jump to the another branch. 

Let us now discuss shortly the experimental situation. 
As we mentioned in Ref^, the BPR exists in the numer- 
ical simulations of the other authors as well. The points 
like breakpoint are presented in the experimental CVC 
of the Refs . 10 ! 11 Recently, new detailed experiment was 
made and the clear manifestation of the BPR structure 
was found.— It's making promising the possibility to ob- 
serve the experimentally the features presented in this 
paper. The following conditions are important for that. 
The measurement should be done with high precision in 
bias current and slow sweeping of the CVC. Because the 
BPR region width is decreasing rapidly with the number 
of junctions 3 -, it should not be too big to observe the BPR 
structure. If the step in the bias current is big enough, it 
makes possible to pass the bifurcation points,— related to 
the parametric resonance in the system. The influence 
of the inductive coupling on the BP and BPR is not in- 
vestigated yet, but if the stack's size is smaller than the 
Josephson penetration depth, we may neglect it. 

In summary, we predict a fine structure in the outer- 
most branch of the current-voltage characteristics in the 
layered superconductors. We establish a correspondence 
between the features of this structure and the character 
of the charge oscillations in the different superconduct- 
ing layers in the stack. We demonstrate that the different 
parts of the BPR correspond to the different dynamics 
of the charge oscillation in the superconducting layers. 
Using an idea about parametric resonance at the break- 
point, we explain the origin of the BPR structure. A 
plasma mode is a stable solution of the system in the BPR 
and this fact might be used in some applications, partic- 
ularly, in high frequency devices such as THz oscillators 
and mixers. We consider as well that the breakpoint phe- 
nomenon might be used to understand the mechanism of 
the powerful coherent radiation from the stack of the IJJ 
in BSCCO. We thank T. Hikihara, I. Kakeya, A. Irie, 
H. J. Lee, K. Kadowaki, N. F. Pedersen, P. Seidel, T. 
Koyama, M. Machida, T. Hatano, N. Plakida for helpful 
discussions. Yu.M.S. thanks Prof. J. Ishikawa and Prof. 
S. Noda for kind hospitality in during his stay at Ky- 
oto university. This research was supported by Russian 
Foundation for Basic Research, grant 08-02-00520-a. 
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